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Flameless  combustion  is  a  well  known  measure  to  reduce  NO*  emissions  in  gas  combustion  but  has  not 
yet  been  fully  adapted  to  pulverised  coal  combustion.  Numerical  predictions  can  provide  detailed  infor¬ 
mation  on  the  combustion  process  thus  playing  a  significant  role  in  understanding  the  basic  mechanisms 
for  pollutant  formation.  In  simulations  of  conventional  pulverised  coal  combustion  the  gasification  by 
C02  or  H20  is  usually  omitted  since  its  overall  contribution  to  char  oxidation  is  negligible  compared  to 
the  oxidation  with  02.  In  flameless  combustion,  however,  due  to  the  strong  recirculation  of  hot  combus¬ 
tion  products,  primarily  C02  and  H20,  and  the  thereby  reduced  concentration  of  02  in  the  reaction  zone 
the  local  partial  pressures  of  C02  and  H20  become  significantly  higher  than  that  for  02.  Therefore,  the 
char  reaction  with  C02  and  H20  is  being  reconsidered.  This  paper  presents  a  numerical  study  on  the 
importance  of  these  reactions  on  pollutant  formation  in  flameless  combustion.  The  numerical  models 
used  have  been  validated  against  experimental  data.  By  varying  the  wall  temperature  and  the  burner 
excess  air  ratio,  different  cases  have  been  investigated  and  the  impact  of  considering  gasification  on 
the  prediction  of  NO  formation  has  been  assessed.  It  was  found  that  within  the  investigated  ranges  of 
these  parameters  the  fraction  of  char  being  gasified  increases  up  to  35%.  This  leads  to  changes  in  the  local 
gas  composition,  primarily  CO  distribution,  which  in  turn  influences  NO  formation  predictions.  Consider¬ 
ing  gasification  the  prediction  of  NO  emission  is  up  to  40%  lower  than  the  predicted  emissions  without 
gasification  reactions  being  taken  into  account. 

©  2009  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

Flameless  combustion  is  of  ever  growing  interest  as  a  measure  to 
reduce  NOx  emissions  from  combustion  processes.  It  is  state  of  the 
art  for  gas  burners  in  the  heat  treatment  industry.  Here,  NOx  emissions 
are  near  the  detection  limit.  In  flameless  combustion  the  reaction 
zone  is  extended  throughout  a  large  volume  of  the  furnace.  Thus, 
flame  formation  with  its  inherent  high  temperature  peaks  is  inhibited 
and  the  temperature  held  and  species  distribution  are  more  homoge¬ 
neous.  This  combustion  mode  is  achieved  by  a  strong  recirculation  of 
flue  gases  to  the  primary  reaction  zone  causing  the  reactants  to  be  sig¬ 
nificantly  diluted.  Recently,  there  is  an  increased  interest  in  applying 
this  technique  to  pulverised  coal  combustion.  It  has  already  been 
shown  experimentally  that  NOx  emissions  of  coal  combustion  could 
be  significantly  reduced  with  such  a  burner  design  [1,2]. 

In  conventional  flames  the  oxygen  concentration  in  the  reaction 
zone  is  high  enough  for  most  of  the  char  to  be  oxidised  by  oxygen. 
This  reaction  is  several  orders  of  magnitude  faster  than  the  reactions 
with  carbon  dioxide  and  water  vapour.  Therefore,  the  gasification 
reactions  are  often  neglected  in  numerical  investigations,  [3-7]. 
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In  flameless  combustion,  however,  the  intense  recirculation  of 
hot  combustion  products  causes  fast  dilution  of  the  02  in  the  reac¬ 
tion  zone.  Thus,  the  partial  pressure  of  C02  and  H20  become  higher 
than  that  of  02  and  therefore  the  oxidation  of  char  is  likely  to  fol¬ 
low  diverse  paths  in  parallel.  Considering  this  in  numerical  simula¬ 
tions  should  lead  to  more  accurate  predictions  of  gas  species 
concentrations.  The  need  for  more  detailed  models  for  char  burn¬ 
out  was  recently  reported  by  Schaffel  et  al.  [8]. 

The  importance  of  consideration  of  the  gasification  reactions 
was  already  observed  for  a  standard  low  NOx  burner  [9].  In  order 
to  identify  the  boundary  conditions  under  which  gasification  reac¬ 
tions  should  be  considered  and  to  get  a  better  understanding  of  the 
pollutant  formation  mechanisms  in  flameless  coal  combustion,  a 
numeric  parameter  study  was  carried  out  to  determine  the  relative 
importance  of  oxygen,  carbon  dioxide  and  water  vapour  as  oxidiser 
for  char  and  the  influence  of  the  changed  temperature  field  and 
species  concentration  on  NO  formation. 

2.  Experimental  setup 

Experiments  were  conducted  in  a  vertical,  cylindrical,  top  fired 
furnace  with  a  diameter  of  400  mm  and  a  total  length  of  4200  mm, 
Fig.  1.  In  this  furnace  the  burner  is  axially  traversable  allowing  for 


0010-2180/$  -  see  front  matter  ©  2009  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved, 
doi:  10.1 01 6/j.combustflame.2009.06.006 


1756 


H.  Stadler  et  al.  /  Combustion  and  Flame  156  (2009)  1755-1763 


0  400 


traversable 
burner 


measurement 

access 


ash  exit 


flue  gas 
exit 


Fig.  1.  Test  facility,  all  dimensions  in  (mm). 


measurements  at  different  distances  from  the  burner  through  a 
single  measurement  plane.  Thus,  for  measurements  of  the  velocity 
the  furnace  is  only  2100  mm  plus  the  actual  distance  of  the  mea¬ 
surement  plane  from  the  burner  whereas  global  measurements 
were  carried  out  with  a  total  furnace  length  of  3900  mm.  During 
in  flame  measurements  burnout  air  was  introduced  700  mm  from 
the  bottom  of  the  furnace,  whereas  during  global  measurements  it 
was  introduced  at  the  measurement  plane  with  the  burner 
1800  mm  above  the  measurement  plane.  Electrical  heating  was 
used  for  regulation  of  the  wall  temperature. 

The  pulverised  coal  burner  employed  in  this  study,  Fig.  2,  was 
developed  according  to  the  design  of  flameless  gas  burners.  The 
predominant  design  criterion  for  a  flameless  combustion  burner 
is  the  avoidance  of  high  temperature  zones  in  order  to  reduce 
the  formation  of  thermal  NO.  In  coal  flames  thermal  NO  makes 
up  only  a  small  fraction  of  the  total  NO  emission  as  fuel  NO  ac¬ 
counts  for  a  major  portion  of  NO  emissions.  However,  experiments 
[2]  have  shown  the  high  NO  reduction  potential  of  flameless  coal 
combustion  which  cannot  be  attributed  to  the  prevention  of  ther¬ 
mal  NO  alone.  In  accordance  with  gas  burners  for  flameless  corn- 


carrier  air  +  coal 


bustion  a  homogeneous  temperature  distribution  is  achieved  by 
good  mixing  of  the  incoming  fresh  mixture  and  the  flue  gas.  A  high 
inlet  momentum  of  the  secondary  air  is  used  to  induce  the  re¬ 
quired  recirculation  of  the  flue  gas  which  is  then  entrained  into 
the  inlet  streams. 

The  dimensions  of  the  investigated  burner  are  given  in  Fig.  2. 
The  coal  is  transported  by  a  carrier  air  flow  and  enters  the  chamber 
through  a  central  inlet.  The  secondary  air  is  injected  through  three 
nozzles  (diameter  4.2  mm)  which  are  arranged  annularly  around 
the  central  coal  jet.  A  sealing  air  stream  was  used  in  the  gap  be¬ 
tween  the  traversable  burner  and  the  furnace  wall  in  order  to  pre¬ 
vent  particles  from  entering  the  volume  above  the  burner.  For  the 
calculation  of  the  burner  excess  air  ratio  this  air  stream  is  also 
considered. 

This  burner  design  was  subject  to  an  experimental  study  inves¬ 
tigating  its  performance  with  different  coals  and  different  stoichi¬ 
ometric  boundary  conditions  [2].  The  arrangement  of  the  nozzles 
allows  recirculated  hot  flue  gases  to  penetrate  between  the  air  jets 
directly  to  the  coal  jet.  Thus,  the  secondary  air  as  well  as  the  coal 
are  diluted  with  flue  gas  before  they  mix.  The  oxygen  concentra¬ 
tion  in  the  reacting  mixture  is  reduced  and  thus  the  peak  flame 
temperature  is  lowered. 

Furthermore,  it  can  be  expected  that  NO  that  is  produced  in  the 
reaction  zone  is  partially  recirculated  together  with  the  flue  gases 
thus  allowing  for  NO  reduction  in  the  burner  vicinity. 

For  the  current  numerical  investigation  reference  data  from 
three  different  experiments  were  considered.  The  boundary  condi¬ 
tions  are  listed  in  Table  1.  The  first  case  at  a  burner  excess  air  ratio 
of  0.7  was  used  to  verify  the  flow  field  as  data  from  velocity  mea¬ 
surements  were  available  for  comparison.  The  other  two  cases 
(burner  excess  air  ratios  0.8  and  0.9)  were  used  for  validation  of 
NO  predictions. 

The  investigated  coal  is  a  blend  from  different  Polish  bitumi¬ 
nous  coals.  Proximate  and  ultimate  analysis  are  given  in  Table  2. 

2.2.  Measurement  technique 

Velocity  measurements  were  carried  out  using  Laser  Doppler 
Anemometry  (LDA).  The  Laser  Doppler  Anemometry  permits 
non-intrusive  measurements  of  axial  and  tangential  velocities  with 
high  spatial  and  temporal  resolution. 

For  measurements  within  the  coal  flame  no  additional  seeding 
of  the  flow  was  necessary  since  the  particle  density  of  coal  and 
ash  particles  is  sufficient  within  the  whole  measurement  volume. 
Due  to  the  small  particle  size  (90%  of  the  coal  mass  is  contained 
in  particles  smaller  than  75  pm;  Dp90  <  75  pm)  it  can  be  assumed 
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Table  1 

Flow  boundary  conditions  and  NO  emissions  determined  from  experiments. 


Burner  EAR 

Coal  flow  00 

Carrier  air  00 

Secondary  air  00 

Sealing  air  00 

Burnout  air  00 

Wall  temperature  (°C) 

Global  NO  emission  (ppm) 

0.7 

1.167  x  10~3 

3.389  x  10~3 

3.699  x  10~3 

2.730  x  10~3 

8.800  x  10~3 

900 

50 

0.8 

1.167  x  10  3 

3.389  x  10~3 

4.561  x  10'3 

2.766  x  10~3 

6.824  x  10~3 

900 

190 

0.9 

1.167  x  10  3 

3.389  x  10~3 

5.783  x  10~3 

2.766  x  10~3 

5.388  x  10”3 

900 

261 

EAR:  excess  air  ratio. 


Table  2 

Proximate  and  ultimate  analysis  in  (mass%)  of  bituminous  coal  blend. 


Proximate  analysis 

Ultimate  analysis 

Water 

0.60 

C 

85.90 

Ash 

3.60 

H 

5.09 

Volatiles 

33.00 

O 

2.94 

Char 

62.80 

N 

1.41 

S 

0.46 

that  the  slip  between  particles  and  gas  is  negligible  so  that  mea¬ 
sured  particle  velocities  directly  reflect  gas  velocity  [10]. 

The  optical  access  at  the  furnace  used  for  the  experiments  re¬ 
quires  a  large  working  distance  between  the  optics  and  the  mea¬ 
surement  position  within  the  flow.  Therefore,  individual  beam 
optics  with  single  beam  adjustment  are  used  as  transmitting  op¬ 
tics.  Hence,  by  increasing  the  beam  distance,  an  angle  of  3.2°  could 
be  retained  so  that  a  reasonable  small  measurement  volume  of 
150  pm  diameter  and  2.5  mm  length  can  be  obtained.  In  order  to 
reduce  the  effort  of  alignment,  a  back  scattering  setup  is  chosen 
where  the  transmitting  and  receiving  optics  are  mounted  on  a  sin¬ 
gle  frame.  The  absolute  values  of  the  individual  particle  velocities 
are  subject  to  a  maximum  error  of  3%  due  to  an  uncertainty  in 
the  determination  of  the  actual  beam  angle  in  the  optical  setup. 

The  velocities  reported  here  are  arithmetic  mean  velocities  of 
the  recorded  individual  particles.  Also  a  95%  confidence  interval 
for  the  calculated  mean  is  given. 

Another  uncertainty  next  to  the  determination  of  the  velocity 
lies  in  the  identification  of  the  exact  measurement  position.  The 
axial  positioning  of  the  burner  occurs  with  an  uncertainty  of 
1  mm.  Also  the  radial  positioning  of  the  optical  system  relative  to 
the  furnace  axis  underlies  an  uncertainty  of  approximately  5  mm. 

3.  Numerical  models 

Numerical  simulations  of  pulverised  coal  combustion  were  car¬ 
ried  out  using  the  CFD  code  Fluent  6.3.26  with  modified  (UDF- 
based)  devolatilisation,  char  oxidation  and  gas  phase  combustion 
submodels.  The  flow  field  was  solved  using  the  k-e-realisable  tur¬ 
bulence  model  with  the  Simple  algorithm  for  velocity-pressure  cou¬ 
pling.  Transport  equations  for  the  mass  fractions  of  seven  different 
gas  species  (volatiles  (CxHyO/SnNm),  water  vapour,  carbon  monox¬ 
ide,  carbon  dioxide,  sulphur  dioxide,  oxygen  and  hydrogen)  were 
calculated  with  source  terms  defined  according  to  the  reaction 
scheme  proposed  below.  Being  the  most  abundant  species  in  the 
gas  mixture,  the  nitrogen  mass  fraction  was  calculated  from  a  mass 
balance.  The  following  kinetic  mechanism  according  to  Toporov 
et  al.  [11]  was  assumed  for  the  combustion  of  volatiles: 

C*HyO,S„Nm  +  g  +  n-£)02  xCO +|h2  +  nS02  +  ™N2  (1) 

and  subsequently: 

CO  +  lo2^C02  (2) 

H2+1q2^H20  (3) 


The  interaction  between  the  turbulence  and  chemical  reactions 
in  the  gas  phase  was  modelled  using  the  Eddy  Dissipation  Concept 
model  according  to  Magnussen  [12],  implemented  and  presented 
by  Heil  et  al.  [13]  with  kinetic  rates  given  in  Table  3. 

Coal  particle  combustion  was  simulated  based  on  a  stochastic 
Lagrangian  procedure  to  track  the  particle  trajectories  in  the  flow 
field.  The  coal  particle  undergoes  decomposition  into  char  and  vol¬ 
atile  material.  The  volatile  material  is  assumed  to  rapidly  form  CO, 
reaction  (1),  from  which  subsequently  C02  is  formed,  reaction  (2), 
whereas  the  char  is  oxidised  in  the  later  stages  of  combustion. 

The  devolatilisation  process  was  modelled  using  a  single  reac¬ 
tion  approach  with  an  activation  energy  of  74  kJ/(mol  I<). 

For  heterogenous  oxidation  of  char  two  different  mechanisms 
have  been  employed.  The  first  mechanism  implies  only  partial  char 
oxidation  by  oxygen: 

Cchar+j02^CO  (4) 

Simulations  with  this  simple  mechanism  were  used  as  compar¬ 
ison  for  the  investigation  of  the  influence  of  the  gasification  reac¬ 
tions.  The  second  mechanism  also  implies  partial  char  oxidation 
by  oxygen  according  to  reaction  (4)  but  extends  further  by  consid¬ 
ering  char  gasification  by  C02  and  H20.  For  modeling  purposes  the 
reaction  of  coal  char  with  C02  and  H20  may  be  simplified  as  the 
reaction  of  carbon  with  C02  and  H20,  respectively: 

Cchar  +  C02  -  2CO  (5) 

Cchar  +  H20  -  CO  +  H2  (6) 

The  Arrhenius  coefficients  used  for  these  reactions  are  given  in 
Table  4.  For  reactions  (4)-(6)  first  order  was  considered. 

However,  one  should  bear  in  mind  that  although  carbon  is  the 
dominant  atomic  species  present  in  coal  char,  its  reactivity  can 
be  quite  different  from  that  of  the  actual  char.  Gasification  rates 
for  a  number  of  coals  and  chars  at  varying  boundary  conditions 
are  available  in  literature  [17-21].  Often  coal  is  more  reactive  than 
pure  carbon  for  a  number  of  reasons.  In  coal  char  various  reactive 
organic  functional  groups  can  enhance  its  reactivity  as  well  as  nat¬ 
urally  occurring  mineral  ingredients  which  enhance  reactivity  cat- 
alytically.  The  actual  reactivity  also  changes  in  time  with 
advancing  burnout  [22].  As  no  experimental  data  are  available 
for  the  coal  investigated  in  this  study,  the  rates  for  graphite  were 
used  which  in  most  cases  are  rather  conservative. 

Rates  given  in  Table  4  are  limited  in  the  temperature  range.  For 
temperatures  out  of  the  given  range  coefficients  valid  for  the  near¬ 
est  regime  were  taken. 


Table  3 

Kinetic  rates  of  gas  reaction,  units  are  m3  kmol  s  J  K. 


Reaction 

A 

E 

Ref. 

(1) 

1.1787  x  107 

5.56  x  107 

Shaw  et  al.  [14] 

(2) 

1.3  x  1011 

1.26  x  108 

Howard  et  al.  [15] 

(3) 

9.87  x  10s 

3.1  x  107 

Hsu  and  Jemcov  [16] 

1758 
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Table  4 

Arrhenius  coefficients  for  heterogeneous  reactions. 


Reaction 

Carbon  type 

MW 

E„/R  (K) 

Temperature  range  (°C) 

Ref. 

(4) 

Char 

5  x  10"3 

8900 

677-1377 

Field  [23] 

(5) 

Graphite 

1.35  x  10  4 

16,300 

850-950 

Smoot  and  Pratt  [24] 

Graphite 

6.35  x  10~3 

19,500 

950-1400 

(6) 

Graphite 

3.19  x  10  1 

25,020 

860-960 

Smoot  and  Pratt  [24] 

Graphite 

1.92  x  10  3 

17,680 

1000-1160 

Pulverised  coal  combustion  was  modelled  with  water  evapora¬ 
tion,  devolatilisation  and  char-burnout  submodels  calculated  in 
parallel. 

The  radiative  heat  source  was  calculated  as  a  function  of  the  lo¬ 
cal  irradiation,  calculated  by  the  Discrete  Ordinates  radiation  mod¬ 
el.  The  local  absorption  coefficient  was  calculated  as  the  sum  of  the 
particle  and  gas  absorption  coefficients.  The  latter  was  based  on 
the  Weighted  Sum  of  Grey  Gases  model  which  considers  the 
absorption  coefficients  of  different  grey  gases  of  the  mixture. 

The  three  dimensional  computational  domain  represents  1/6  of 
the  furnace  volume  with  symmetric  boundary  conditions.  The 
numerical  grid  contains  370,000  cells.  The  emissivities  of  the  furnace 
and  burner  walls  were  taken  to  be  s  =  0.7  and  s  =  0.3,  respectively. 

In  flow  conditions  were  determined  with  a  preceding  simula¬ 
tion  of  the  burner  pipe  flow.  Thus,  boundary  conditions  for  velocity 
and  temperature  distribution  across  the  inlet  were  obtained. 

4.  Validation  of  simulations  with  experiments 

Data  from  LDA  velocity  measurements  were  used  to  verify  the  in¬ 
let  boundary  conditions  of  the  numerical  simulation.  Measurements 
had  been  carried  out  at  several  planes  25  mm,  50  mm,  75  mm  and 
1 00  mm  from  the  burner.  On  these  planes  radial  profiles  were  taken. 
The  line  of  measurements  was  inclined  to  the  axis  from  the  primary 
air  to  the  secondary  air  jet  by  about  16°.  Thus,  no  data  from  within 
the  secondary  air  jet  in  burner  vicinity  is  available. 

In  Fig.  3  experimental  and  numerical  data  are  compared.  At  a 
distance  of  25  mm  from  the  burner  (uppermost  graph)  the  exper¬ 
imental  data  reflects  the  central  carrier  air/coal  jet  with  a  maxi¬ 
mum  velocity  of  around  21  m/s.  The  secondary  air  does  not 
extend  to  the  line  of  measurements  at  this  plane.  Yet,  a  minor 
influence  of  the  secondary  air  jet  can  be  seen  in  the  simulation. 
This  difference  might  be  due  to  the  uncertainty  in  axial  positioning 
of  the  burner  described  above. 

At  a  distance  of  50  mm  (second  from  top  in  Fig.  3)  the  secondary 
air  jet  has  spread  far  enough  to  be  detected  along  the  line  of  mea¬ 
surements.  It  appears  as  a  second  peak  at  a  radial  distance  of  20 
mm  from  the  axis.  Its  maximum  velocity  is  about  16  m/s.  The  mea¬ 
sured  minimum  between  the  peaks  is  not  as  low  as  in  the  numerical 
simulation  which  is  due  to  a  lower  resolution  of  the  experiments 
and  the  length  of  the  measurement  volume  of  about  2.5  mm.  Data 
from  within  this  measurement  volume  is  integrated  so  that  struc¬ 
tures  smaller  than  its  length  cannot  be  fully  resolved.  Qualitatively, 
data  from  the  planes  at  75  and  100  mm  (lower  two  graphs  in  Fig.  3) 
look  the  same,  only  peak  values  differ  slightly.  The  maximum  veloc¬ 
ities  are  underpredicted  by  about  3  m/s.  The  value  of  the  local  min¬ 
imum  observed  in  between  the  two  jets  increases  with  increasing 
axial  distance  as  the  two  jets  merge.  Experimental  and  numerical 
data  agree  well  except  for  an  offset  of  the  experimental  data  of 
approximately  5  mm.  This  offset  is  due  to  a  minor  misalignment 
of  the  optical  system. 

The  agreement  obtained  between  experimental  and  numerical 
data  shows  that  the  chosen  k-e-realisable  turbulence  model  and 
the  related  inlet  boundary  conditions  are  appropriate  for  the  inves¬ 
tigated  flameless  pulverised  coal  combustion.  Therefore,  this  case 


Fig.  3.  Comparison  of  experimental  and  numerical  data  at  distances  of  25-100  mm 
from  the  burner. 

was  used  as  a  basis  for  a  case-study  based  investigation  of  the 
influence  of  gasification  reactions  on  NO  emission  predictions. 
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4A.  NO x-models 

Standard  Fluent  NO-postprocessing  models  considering  thermal 
NO,  prompt  NO,  fuel  NO  formation  and  NO  reburning  with  volatiles 
were  used.  Thermal  NO  was  calculated  according  to  the  extended 
Zeldovich  mechanism  with  a  partial  equilibrium  approach  for  the 
radicals  O  and  OH.  For  prompt  NO  volatiles  were  considered  the 
fuel  species. 

One  third  of  the  total  fuel  nitrogen  is  assumed  to  be  released 
with  the  volatiles  and  the  remaining  two  thirds  are  released  to¬ 
gether  with  the  char  burnout  since  the  volatiles  make  up  about 
one  third  of  the  coal  mass  and  the  char  makes  up  two  thirds  of 
the  coal  mass.  For  both  pathways  the  same  proportioning  of  the 
intermediates  was  chosen.  It  was  assumed  that  60%  of  the  nitrogen 
are  released  via  the  intermediate  HCN  and  10%  via  the  intermedi¬ 
ate  NH3.  The  remaining  nitrogen  is  released  directly  as  NO.  Reac¬ 
tions  considered  for  NO  formation  and  those  leading  to  the 
reduction  of  NO  are  the  following: 


HCN  +  02 - 

+  NO+-- 

(7) 

HCN  +  NO  - 

-N2  +  ..- 

(8) 

NH3  +  02  - 

>  NO  +  •  •  ■ 

(9) 

NH3  +  NO  - 

-  N2  +  •  •  • 

(10) 

The  N20  path  has  been  completely  omitted  since  preliminary 
simulations  showed  its  influence  to  be  insignificant  throughout 
the  range  of  parameters  relevant  to  this  study. 

In  the  temperature  range  between  1300  °C  and  1800  °C  NO 
reburning  with  volatiles  was  incorporated.  A  partial  equilibrium 
approach  was  considered  for  the  reburning  species.  The  volatiles 
have  a  C/H  ratio  of  approximately  2.6.  Therefore,  CH3  was  consid¬ 
ered  the  equivalent  reburn  fuel. 

These  parameters  were  validated  against  global  NO  emissions 
measured  during  two  different  experiments  at  a  wall  temperature 
of  900  °C  and  a  burner  excess  air  ratio  (EAR)  of  0.8  and  0.9,  Table  5. 
Predictions  at  a  burner  excess  air  ratio  of  0.8  correspond  very  well 
with  experimental  data,  at  a  burner  excess  air  ratio  of  0.9  NO  emis¬ 
sions  are  slightly  underpredicted. 

5.  Results  and  discussion 

A  number  of  simulations  has  been  carried  out  with  the  wall 
temperature  set  between  900  °C  and  1500  °C  and  at  burner  excess 
air  ratios  between  0.7  and  0.9  as  listed  in  Table  6. 

The  key  features  of  flameless  combustion  are  well  resolved.  The 
intense  recirculation  causes  the  02  concentration  in  the  reaction 
zone  to  be  lower  than  the  C02  concentration.  This  high  dilution 
is  considered  favourable  for  low  NOx  combustion.  Fig.  4  shows 
the  predicted  temperature  distribution  with  and  without  gasifica¬ 
tion  reactions  at  a  burner  excess  air  ratio  of  0.9  and  a  wall  temper¬ 
ature  of  1500  °C.  The  temperature  throughout  the  whole 
combustion  zone  is  characterised  by  small  gradients  and  low  peak 
temperatures  compared  to  the  wall  or  outlet  temperature.  The  pre¬ 
dictions  show  that  the  gas  temperature  distribution  is  independent 
of  the  chosen  char  oxidation  model  and  it  is  homogeneous  in  the 
primary  combustion  zone.  The  difference  in  the  obtained  maxi¬ 
mum  temperatures  between  the  different  cases  does  not  exceed 
20°. 


Table  5 

NO  emission  in  experiments  compared  to  numerical  predictions. 


Table  6 

Boundary  conditions  investigated. 


Burner  EAR  Wall  temperature  (°C) 
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Fig.  4.  Temperature  distribution  in  the  furnace  as  simulated  without  gasification 
reactions  (left)  and  with  gasification  reactions  (right)  at  a  burner  excess  air  ratio  of 
0.9  and  a  wall  temperature  of  1500  °C. 
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However,  within  this  study  it  was  found  that  the  char  gasification 
reactions  can  have  significant  influence  on  overall  char  oxidation. 
Figs.  5  and  6  show  the  share  of  CO  production  from  the  three  reac¬ 
tions  considered  for  char  oxidation  at  a  burner  excess  air  ratio  of 
0.7  (Fig.  5)  and  0.9  (Fig.  6),  predicted  at  different  wall  temperatures. 
An  increase  in  temperature  increases  the  fraction  of  char  being  gas¬ 
ified.  At  a  wall  temperature  of  1 500  °C  and  a  burner  excess  air  ratio  of 
0.7  up  to  30%  of  the  char  are  gasified  by  C02  according  to  reaction  (5) 
and  up  to  5%  are  gasified  by  H20,  reaction  (6).  At  a  burner  excess  air 
ratio  of  0.9  up  to  20%  are  gasified  by  C02  and  3%  by  H20.  At  a  given 
wall  temperature  the  fraction  of  char  being  gasified  by  C02  or  H20 
at  a  burner  excess  air  ratio  of  0.7  is  always  higher  than  at  a  burner  ex¬ 
cess  air  ratio  of  0.9.  It  shows  that  the  02  concentration  is  the  second 
parameter  influencing  the  rates  of  gasification.  At  lower  local  02  con¬ 
centrations  a  bigger  part  of  the  char  is  gasified. 

For  the  given  reaction  rates  for  graphite  (Table  4)  the  reaction  of 
C02  always  yields  about  5-10  times  as  much  CO  as  the  reaction 
with  H20. 

The  gasification  reactions  can  change  the  species  concentrations 
significantly.  With  gasification  reactions  a  bigger  portion  of  the 
char  is  oxidised  resulting  in  higher  CO  concentrations  in  the  pri¬ 
mary  combustion  zone  before  staging  air  injection.  Since  the  CO 
concentration  is  increased,  oxygen  is  more  rapidly  being  consumed 
by  means  of  reaction  (2).  This  in  turn  reduces  the  oxygen  concen¬ 
tration  in  the  substoichiometric  zone  further.  At  the  lowest  inves- 
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Fig.  5.  Relative  importance  of  gasification  reactions  as  CO  source  at  burner  excess 
air  ratio  of  0.7. 
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Fig.  6.  Relative  importance  of  gasification  reactions  as  CO  source  at  burner  excess 
air  ratio  of  0.9. 


tigated  temperatures  (900  °C)  and  at  a  burner  excess  air  ratio  of  0.9 
the  maximum  CO  concentration  obtained  in  simulations  with  gas¬ 
ification  reactions  is  about  10%  higher  than  in  the  case  without. 
Fig.  7  shows  the  CO  distribution  at  a  wall  temperature  of  1500  °C 
on  the  left  side  as  simulated  without  gasification  reactions  and 
on  the  right  side  with  consideration  of  the  gasification  reactions. 
Here,  the  obtained  maximum  CO  concentration  in  simulations  with 


Fig.  7.  Distribution  of  mass  fraction  of  CO  in  the  furnace  as  simulated  without 
gasification  reactions  (left)  and  with  gasification  reactions  (right)  at  a  burner  excess 
air  ratio  of  0.9  and  a  wall  temperature  of  1500  °C. 
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gasification  reactions  is  about  twice  as  high  as  in  the  case  without. 
The  volume  with  high  CO  concentration  is  enlarged  from  a  local¬ 
ised  zone  near  the  burner  to  the  whole  primary  combustion  zone 
up  to  the  injection  of  the  staging  air.  In  the  simulation  without  gas¬ 
ification  reactions  (left)  high  values  show  up  only  in  the  burner 
vicinity  where  the  volatiles  are  burned.  The  CO  distribution  as  sim¬ 
ulated  with  gasification  reactions  (on  the  right)  features  this  same 
peak  but,  as  all  oxygen  is  consumed  for  the  oxidation  of  CO,  the  CO 
concentration  starts  increasing  up  to  the  inlet  of  the  staging  air. 
The  path  of  the  recirculation  also  becomes  obvious  as  a  zone  of 
high  CO  concentration.  It  extends  along  the  wall  up  to  the  inlet 
of  the  sealing  air  where  the  CO  is  oxidised.  Thus,  a  major  portion 
of  the  recirculation  occurs  in  a  highly  reducing  atmosphere  favour¬ 
ing  low  NO  production. 

This  in  turn  has  an  influence  on  the  prediction  of  pollutant  for¬ 
mation.  As  the  oxygen  concentration  decreases,  the  concentrations 
of  the  NO-intermediate  products  HCN  and  NH3  increase  because 
reactions  (7)  and  (9)  are  impeded.  The  oxygen  is  more  rapidly  con¬ 
sumed  for  CO  oxidation  (2)  as  for  the  oxidation  of  the  intermedi¬ 
ates.  The  increased  intermediate  concentration  in  turn  leads  to 
higher  reaction  rates  for  the  reduction  routes,  reactions  (8)  and 
(10)  resulting  in  a  lower  NO  concentration. 

An  accurate  prediction  of  the  CO  distribution  becomes  even 
more  important  when  its  influence  on  the  radical  pool  is  consid¬ 
ered.  According  to  Brouwer  et  al.  [25]  and  Kristensen  et  al.  [26]  un¬ 
der  selective  non  catalytic  reduction  (SNCR)  conditions  a  high  CO 
concentration  is  beneficial  for  the  reduction  of  NO  because  it  en¬ 
larges  the  radical  pool  by  its  reaction  with  OH  yielding  C02  and 
H  and  subsequent  chainbranching  reactions  of  H.  Under  conditions 
of  flameless  coal  combustion  a  similar  influence  of  high  CO  concen¬ 
tration  might  exist  but  further  work  is  needed  for  clarification. 
With  the  global  models  used  in  this  study  this  effect  was  omitted. 

5.2.  Results  from  NO  simulations 

Fig.  8a  shows  the  calculated  global  NO  emissions  at  different 
wall  temperatures  for  both  cases  -  with  and  without  gasification 
reactions.  For  a  burner  excess  air  ratio  of  0.9  NO  emission  drops 
continuously  with  increasing  wall  temperature  up  to  a  wall  tem¬ 
perature  of  1300  °C.  At  higher  temperatures  a  strong  increase  in 
NO  emissions  can  be  seen.  For  a  burner  excess  air  ratio  of  0.7  the 
trend  is  not  so  obvious.  Between  900  °C  and  1100  °C  the  NO  emis¬ 
sion  increases  and  drops  again  towards  1400  °C.  A  further  temper¬ 
ature  increase  leads  to  an  increased  NO  emission  again.  This 
behaviour  is  the  sum  of  the  individual  mechanisms  and  can  be  ex¬ 
plained  when  the  origin  of  NO  is  considered. 

Fig.  8b  shows  the  NO  produced  via  the  fuel-NO  path.  The  NO 
models  used  for  the  current  analysis  capture  the  trend  of  lower 
fuel-NO  production  at  higher  temperatures  as  it  was  reported  in 
several  experimental  studies  [2,27,28].  With  higher  temperatures 
less  fuel-NO  is  produced  because  the  reaction  of  the  intermediates 
HCN  and  NH3  with  NO  to  form  N2,  reactions  (8)  and  (10),  are  faster 
than  the  oxidation  reactions  (7)  and  (9)  especially  since  the  fuel  is 
released  in  a  highly  substoichiometric  region.  This  difference  be¬ 
comes  more  significant  with  increased  temperature. 

At  temperatures  above  1300  °C  the  predicted  fuel-NO  emission 
increases.  This  effect  is  due  to  the  exponential  increase  of  reaction 
rates  with  temperature  but  does  not  necessarily  reflect  real 
behaviour  as  the  rate  coefficients  used  do  not  cover  this  high  tem¬ 
perature  range.  Similar  predictions  and  comparison  with  experi¬ 
mental  data  have  been  obtained  by  Alzueta  et  al.  [29]. 

Fig.  8c  shows  the  contribution  of  the  reburning  reactions.  The 
magnitude  of  the  reburning  rates  is  of  the  same  order  as  the 
fuel-NO  production  compensating  for  most  of  the  production. 
However,  with  increasing  temperature  the  amount  of  NO  that  is  re¬ 
burned  is  reduced  and  seems  to  stagnate  between  1200  °C  and 
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A  EAR  0.9  without 


•  EAR  0.7  with 
■  EAR  0.8  with 
A  EAR  0.9  with 


wall  temperature  [°C] 


Fig.  8.  Total  predicted  NO  emission  and  distribution  among  mechanisms  of  NO 
production  and  reduction. 


1500  °C  causing  the  overall  emission  of  NO  to  rise  with  tempera¬ 
ture.  A  possible  reason  for  the  stagnation  is  that  the  volatiles, 
which  are  considered  the  reburning  species  in  the  current  investi¬ 
gation,  are  consumed  faster  at  higher  temperatures  so  that  the 
slower  reaction  with  NO  is  impeded. 
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Fig.  9.  Relative  difference  of  NO  emission  with  and  without  gasification  reactions 
NQno~N°w  (NOwo:  predicted  NO  emission  without  gasification  reactions,  NOw: 
predicted  NO  emission  with  gasification  reactions)  over  share  of  oxidation  by  C02 
gasification  reaction. 


Fig.  8d  shows  the  NO  production  via  the  thermal  path.  Up  to 
1400  °C  the  overall  production  is  negligible  compared  to  fuel-NO 
production  (less  than  5%).  This  is  a  direct  influence  of  the  flameless 
combustion  mode  where  peak  temperatures  are  kept  low  and  thus 
thermal  NO  formation  is  reduced.  At  a  wall  temperature  of  1500  °C 
a  considerable  amount  of  NO  is  produced  via  the  thermal  NO 
mechanism,  however  NO  formation  via  the  fuel-NO  path  yields 
about  five  times  as  much  NO.  In  experiments  in  a  similar  furnace 
[2]  but  with  a  different  coal  (lignite  instead  of  bituminous  coal) 
it  was  found  that  about  one  third  of  total  NO  emissions  in  flameless 
combustion  are  produced  via  the  thermal  NO  mechanism  at  a  wall 
temperature  of  1300  °C. 

Mancini  et  al.  [30]  reported  that  the  prompt  NO  route  is  domi¬ 
nant  in  the  fuel  rich  primary  reaction  zone  of  flameless  combustion 
of  gaseous  fuels,  however,  the  NO  mechanism  used  in  the  current 
study  showed  only  a  negligible  contribution  throughout  the  whole 
temperature  range  investigated. 

Calculated  NO  values  without  gasification  reactions  are  always 
higher  than  in  the  case  when  the  gasification  reactions  are  consid¬ 
ered.  This  is  in  accordance  with  Yamamoto  et  al.  [9].  The  difference 
increases  with  temperature  and  decreases  with  the  burner  excess 
air  ratio.  In  Fig.  9  the  difference  in  NO  prediction  is  plotted  against 
the  share  of  C02  gasification  reaction  in  total  char  oxidation.  This 
figure  points  up  the  influence  of  the  gasification  reactions  on  NO 
prediction.  Up  to  5%  oxidation  of  char  by  C02  the  difference  of 
the  NO  prediction  between  the  simulation  with  and  without  gasi¬ 
fication  reaction  is  negligible.  In  this  case  the  species  concentra¬ 
tions  do  not  differ  significantly.  At  an  increased  gasification 
fraction  of  7.5%  there  is  already  a  10%  difference  between  the  pre¬ 
dictions.  This  difference  increases  to  about  40%  at  a  level  of  around 
45%  gasification  of  the  char  by  C02. 

6.  Conclusions 

A  numerical  investigation  on  the  influence  of  char  gasification 
reactions  on  NO  emission  predictions  has  been  made.  The  numer¬ 
ical  model  has  been  validated  against  in-flame  measurements  of 
the  velocity  and  global  NO  emission  measurements.  A  case  study 
approach  was  used  for  the  assessment  of  the  relative  importance 
of  gasification  reactions.  The  wall  temperatures  assumed  in  this 
study  were  extended  beyond  those  values  applicable  to  standard 
combustion  systems  in  order  to  capture  the  trend.  Flowever,  the 
char  reactivity  strongly  depends  on  the  coal  type  and  results  ob¬ 
served  at  high  temperatures  in  this  study  are  likely  to  be  observa¬ 


ble  at  lower  temperatures  for  coals  which  are  more  reactive  than 
the  graphite  underlying  the  current  study. 

The  use  of  kinetic  data  for  carbon  burnout  that  is  derived  for  the 
coal  investigated  in  numerical  simulations  will  lead  to  quantitative 
improvements  of  the  predictions. 

This  work  has  shown  that  in  flameless  pulverised  coal  combus¬ 
tion  a  non  negligible  portion  of  char  can  be  oxidised  by  C02  and 
FI20.  Considering  these  reactions  in  simulations  of  flameless  com¬ 
bustion  has  negligible  influence  on  the  predicted  temperature  dis¬ 
tribution  but  results  in  different  species  concentration  fields. 
Hence,  the  prediction  of  pollutant  formation  can  alter  significantly. 
In  particular  the  findings  are  as  following: 

1.  The  fraction  of  char  oxidised  by  C02  and  H20  increases  with 
temperature  and  decreases  with  higher  burner  excess  air  ratio. 

2.  Gasification  reactions  increase  the  burnout  in  the  primary  com¬ 
bustion  zone. 

3.  C02  and  H20  gasification  lead  to  an  increased  CO  concentration 
in  the  primary  combustion  zone  thus  further  reducing  the  02 
concentration. 

4.  The  reduced  oxygen  concentration  leads  to  a  decrease  of  the 
concentration  of  NO  intermediate  species  in  the  primary  reac¬ 
tion  zone. 

5.  Less  oxygen  and  higher  NO-intermediate  concentrations  result 
in  higher  reduction  of  NO. 

Based  on  these  observations  it  can  be  concluded  that  the  design 
of  a  flameless  coal  combustion  burner  should  promote  intense 
mixing  of  recirculated  hot  combustion  products  and  incoming 
fresh  reactants  thus  allowing  simultaneous  fast  02  dilution  and 
fuel  heat-up  which  as  a  result  favours  char  gasification  with  the 
benefit  of  enhanced  NO  reduction.  This  design  also  leads  to  a 
homogeneous  temperature  distribution  which  is  considered 
advantageous  for  reduced  thermal  NO  formation. 

Although  the  NO  model  used  in  this  study  was  able  to  capture 
the  correct  trend  and  provide  reasonable  predictions,  a  more 
sophisticated  NO  model  could  give  more  accurate  predictions. 
However,  the  most  recent  NO  models  available  in  literature  have 
not  yet  been  tested  and  validated  for  flameless  coal  combustion. 
Thus,  the  findings  from  this  work  can  be  used  as  a  basis  for  valida¬ 
tion  of  newer  NO  models  for  flameless  coal  combustion. 
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